Theory of hopping conduction in arrays of doped semiconductor nanocrystals 
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The resistivity of a dense crystalline array of semiconductor nanocrystals (NCs) depends in a 
sensitive way on the level of doping as well as on the NC size and spacing. The choice of these 
parameters determines whether electron conduction through the array will be characterized by 
activated nearest-neighbor hopping or variable-range hopping (VRH). Thus far, no general theory 
exists to explain how these different behaviors arise at different doping levels and for different types 
of NCs. In this paper we examine a simple theoretical model of an array of doped semiconductor 
NCs that can explain the transition from activated transport to VRH. We show that in sufficiently 
small NCs, the fluctuations in donor number from one NC to another provide sufficient disorder 
to produce charging of some NCs, as electrons are driven to vacate higher shells of the quantum 
confinement energy spectrum. This confinement-driven charging produces a disordered Coulomb 
landscape throughout the array and leads to VRH at low temperature. We use a simple computer 
simulation to identify different regimes of conduction in the space of temperature, doping level, and 
NC diameter. We also discuss the implications of our results for large NCs with external impurity 
charges and for NCs that are gated electrochemically. 
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I. INTRODUCTION 

Arrays of semiconductor nanocrystals (NCs) have 
great promise for optoelectronic and photovoltaic de- 
vices. The usefulness of NC arrays comes from the ability 
to tune both their optical properties - generally by choos- 
ing the size or shape of NCs [ll l2| - and their electronic 
properties - usually through the addition of dopants or 
surface ligands that control the spacing between NCs 
[3,13 . Recent experiments have demonstrated that dense, 
crystalline arrays of spherical semiconductor NCs can be 
reliably produced with diameter in the range 4-10 nm 
and with less than 5% dispersion [HQ- Thus, optoelec- 
tronic or photovoltaic devices made from NCs can be 
designed to operate precisely in any chosen region of the 
optical spectrum. 

From a practical standpoint, however, the development 
of NC-based devices is slowed by the high resistivity of 
the NC arrays. In their undoped state, semiconductor 
NCs are insulators, and in order to reduce their large re- 
sistivity it is necessary to bring additional electrons (or 
holes) to the NCs either^through chemical doping 6] or 
electrochemical gating [7|. In this article we focus pri- 
marily on the former, although we comment on electro- 
chemical gating at the end of the paper. 

In particular, we consider the case where each NC is 
made from a semiconductor that is heavily-doped, for 
example, by donor impurities. In this case all donor elec- 
trons reside in the conduction band of the NC. In order 
to conduct across the array, these electrons must tunnel 
between NCs under the high barrier associated with the 
insulator (such as the ligands shown in Fig. [T]) that fills 
the space between them. 

In the presence of even a relatively small amount of 
disorder in the array, the large tunneling barriers im- 
ply that donor electrons experience Anderson localiza- 
tion due to fluctuations in the electron energy from one 
NC to another Q. In this situation conduction pro- 




FIG. 1. (Color online) Schematic drawing of spherical semi- 
conductor NCs (large, light-colored circles) with diameter D 
arranged in a crystalline lattice with lattice constant D' . Each 
NC is coated in a thin layer of insulating ligands (curvy lines) 
that maintain a separation d = D' — D between NCs and pre- 
vent them from sintering. Each NC has a random number of 
donors in its interior (small, black circles). 



ceeds only by phonon-assisted tunneling, or "hopping", 
between localized electron states. This hopping is a 
thermally-activated process in which electron tunneling 
occurs simultaneously with the absorption or emission 
of a phonon whose energy accounts for the difference be- 
tween the initial and final electron states. (While metallic 
conduction through the array is in principle possible, and 
has been reported [9| , it requires the characteristic disor- 
der energy in the system to be smaller than the hopping 
integral t between neighboring NCs. Since t decays ex- 
ponentially with the separation d between NCs and with 
the height of the tunneling barrier them, the condition 



for metallic conductivity is difficult to meet, and in this 
paper we assume that electron conduction proceeds by 
hopping.) 

If one assumes that in the global ground state of the ar- 
ray all NCs are neutral, then hopping transport requires 
an electron to be thermally excited to jump from one 
neutral NC to another. This process produces two op- 
positely charged NCs, each of which has a corresponding 
Coulomb self-energy £c — e^ / k,D, where k, is the effective 
dielectric constant of the NC array and D is the NC diam- 
eter. This charging energy plays the role of an activation 
energy for resistivity in the case where all NCs are neu- 
tral in the global ground state. Equivalently, one can say 
that the distribution of electron ground state energies, 
or the "density of ground states" (DOGS) of NCs, has 
a gap of width 2ec centered at the electron Fermi level. 
As a result, the resistivity p follows the Arrhenius law: 
Inp oc Ec/kBT, where fcsT is the thermal energy. We 
emphasize that the activation energy for hopping con- 
duction is sensitive only to the ground state energies of 
electrons and holes that are added to NCs. For this rea- 
son when calculating the resistivity it is sufficient to con- 
sider the DOGS, which does not include excited electron 
states with additional kinetic energy. 

In experiments, however, one often observes a temper- 
ature dependence of the resistivity that is different from 
simple activation: Inp ex T"'"', with the temperature ex- 
ponent 7 < 1. Such "stretched exponential" behavior is 
believed to be possible only if the disorder is so strong 
that a substantial fraction of NCs is charged in the global 
ground state. Such charging creates a random Coulomb 
potential landscape that shifts up and down the electron 
energy spectra at different NCs. As a result of this shift- 
ing, the gap in the DOGS is smeared and filled. This 
smearing means that some electron states have energies 
very close to the Fermi level, and as a result one can find 
a pair of empty and filled electron states separated by 
an energy As that is much smaller than Ec- At small 
temperature ksT ^ Ec, it is hopping between such pairs 
that are close in energy that dominates the conduction. 

Of course, for small Ae the typical separation r be- 
tween the corresponding NC pair is much larger than 
the spacing D' between neighboring NCs. Thus, at 
small temperature T electron conduction relies on tun- 
neling between distant NCs. To understand how such 
long-range tunneling is possible, consider first the tun- 
neling of an electron between nearest-neighboring NCs. 
When the electron tunnels through the insulating gap 
of thickness d between NCs, it accumulates an action 
hd/ a, where a is the decay length of the electron wave- 
function outside of the NC. Thus, the tunneling ampli- 
tude between nearest neighbors is suppressed by a factor 
~ exp[— d/a]. On the other hand, when an electron tun- 
nels to a NC at a distance x ^ D' , the path of least action 
for the electron is to travel primarily through nearest- 
neighboring NCs, making hops only through the small 
gaps between neighbors and thereby accumulating an ac- 
tion ~ h{d/a){x/ D'), plus an additional much smaller 



term corresponding to action accumulated across the in- 
terior of each NC. Thus, the tunneling amplitude to the 
distance x is suppressed by a factor ^ exp[— a;(i/D'a]. 
The exponential decay of the tunneling amplitude is de- 
scribed by defining the localization length ^, such that 
tunneling between NCs with separation r is suppressed 
by the factor exp[— 2r/^]. By the argument above, one 
cannot simply equate ^ with a, but rather ^ ^ aD' / d 3> a 
|1(]| . It is this enhanced localization length, made possi- 
ble by tunneling through intermediate NCs, that allows 
for long-range hopping. In the remainder of this paper, 
we consider the limit where d and a are both very small 
compared to the NC diameter, so that D' ^ D while f 
remains finite. 

If the temperature T is made increasingly small, the 
corresponding energy difference Ae of electron hops be- 
comes increasingly small due to the scarcity of avail- 
able high-energy phonons, and as a result the typical 
hop length increases. Such behavior is known as vari- 
able range hopping (VRH), and is responsible for the 
stretched exponential behavior 7 < 1 in the resistivity. 
When the DOGS is constant near the Fermi level, the re- 
sistivity follows the Mott law of VRH ll]: Inp ex T''^/^. 
However, in systems where the long-ranged Coulomb po- 
tential is not screened, electron correlation effects pro- 
duce a DOGS that vanishes quadratically with energy 
at the Fermi level [IjI- Such a vanishing DOGS results 
in the Efros-Shklovskn (ES) law of VRH: Inp oc T^^/^. 
In principle, all three of these conduction behaviors — 
Arrhenius (7 = 1), Mott VRH (7 = 1/4), and ES VRH 
(7 — 1/2) — are possible in arrays of semiconductor NCs, 
depending on the magnitude and type of disorder present. 
In this paper we focus our description on the fundamen- 
tal role played by inherent fiuctuations in donor number 
among doped NCs. 

Experiments probing the resistivity of NC arrays have 
reported that the resistivity depends in a sensitive and 
qualitative way on the level of doping [7|. Specifically, 
as the average number v of dopant electrons per NC is 
varied, the dependence of the resistivity p on the temper- 
ature T changes between Arrhenius-type activated con- 
duction (7 = 1) and VRH (7 < 1). VRH has been re- 
ported in a variety of granular semiconductor systems 
0, 0, [13, [iJI , but thus far there is no general theory to 
explain how these different types of conduction can co- 
exist and why they appear in particular ranges of the 
electron "filling factor" v. 

In this paper we present such a theory, based on a first- 
principles description of the ground state arrangement of 
electrons within an array of doped NCs. We focus on 
a simple model of identical spherical NCs that are cov- 
ered by a thin layer of insulating ligand (or some other 
insulator) and arranged in an ideal crystalline lattice, as 
depicted in Fig. [T] We show that the presence of fluctua- 
tions in donor number between different NCs is sufficient 
to produce charging of NCs, which results in a disordered 
Coulomb landscape that encourages VRH. This charging 
is driven by the large gaps between shells of the elec- 



tron quantum energy spectrum in NCs with large Bohr 
radius as- Specifically, these inter-shell gaps drive elec- 
trons to depart from NCs with a large number of donors, 
where maintaining electroneutrality would require plac- 
ing electrons in higher quantum energy shells, and reside 
instead on nearby NCs with small donor number. In this 
way some NCs spontaneously acquire a positive or neg- 
ative charge, and it is this charging that leads to VRH 
when the temperature is not too large. 

Using this model, we explain how the different regimes 
of resistivity observed in experiment arise based on the 
interplay between the charging spectrum of NCs, the 
long-ranged Coulomb interactions between charged NCs, 
and the discrete quantum energy levels of confined elec- 
trons. We supplement our theory with a simple computer 
simulation, which we use to calculate the DOGS and the 
resistivity. 

Our main result is that VRH appears when the average 
number i^ of electrons per NC, the NC diameter D, and 
the temperature T satisfy the following three conditions: 

(i) 1^ Z 0-6, 



5 - □ □[ 



(ii) D < 34KaB/KNCi and 
(iii) keT < 0.5e2^/KL>2. 



Here, knc is the internal dielectric constant of NCs. 
When these three conditions are satisfied, the resistivity 
follows the ES law. In situations where any of the three 
criteria is not met, the conduction is activated. This re- 
sult is depicted at low temperature, ksT ^ e^^/ nD^ , in 
the phase diagram of Fig. [2] 

The remainder of this paper is organized as follows. 
In Sec. In] we define the theoretical model to be studied. 
Sec, mil describes our computer simulation, including our 
methods for numerically calculating the DOGS and re- 
sistivity. Results are presented in Sec. IIVI along with a 
discussion of why Arrhenius and VRH resistivity appear 
in particular regimes oi v, £), and T . We also discuss 
interesting features of the DOGS in this model, includ- 
ing the appearance of "refiected Coulomb gaps" at either 
side of the Fermi level. 

In Sec. |V] we discuss the implications of our theory for 
large NCs with external impurity charges. We present a 
modified model appropriate for this case and we arrive 
at a single condition for VRH associated with the con- 
centration of external impurity charges. Sec. IVII presents 
some speculation on how our results can be applied to 
electrochemical gating of NC arrays using ionic liquids, 
and this is followed by concluding remarks in Sec. IVIII 



II. MODEL OF NC ARRAYS WITH RANDOM 
NUMBER OF DOPANTS 

In this paper our goal is to describe the resistivity of 
a dense array of semiconductor NCs and capture its de- 
pendence on doping level, temperature, and NC diame- 
ter. To this end we adopt the following simplified theo- 
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FIG. 2. (Color online) Phase diagram indicating regimes of 
activated and ES resistivity as a function of doping level u and 
the dimensionless quantum energy gap A = 2Q.64:KaB / k^cD 
at low temperature ksT <C e^^jnD^ . Symbols correspond 
to simulated systems: filled (light blue) circles indicate sys- 
tems that exhibited ES resistivity and open squares indicate 
systems that exhibited activated resistivity. The simulation 
method is described in detail in Sec. lIIII The thick (red) curve 
is an approximate boundary between these two regimes, which 
are labeled "ES" and "A", respectively. Dashed, horizontal 
lines indicate the value of A corresponding to Si NCs with 
D = 5 nm (as in Ref. 01) and to CdSe NCs with D = 6.2 
nm (as in Ref. -7]). This phase diagram is discussed more 
thoroughly in Sec. IIVI 



retical model. We consider NCs to be identical spheres 
of diameter D with large internal dielectric constant 
knc ^ 'S- These spheres are arranged in a regular, three- 
dimensional (3D) lattice, with each lattice site i located 
at the center of a NC. For simplicity, we consider a cubic 
lattice with lattice constant D' just barely larger than D, 
so that d '^ D (see Fig. [T|). Our choice of a cubic lattice 
does not qualitatively affect any of the results we present 
below. 

We further assume that the radius D/2 of the NCs is 
comparable to or smaller than the effective electron Bohr 
radius as = W'K^drae^ of the semiconductor, where e is 
the electron charge and m is the effective electron mass. 
As an example, NCs made from Si have as ~ 2.4 nm; for 
CdSe NCs, a_B ~ 5 nm. Under this condition the wave- 
function of a donor electron is extended across the entire 
volume of a NC, rather than localized around a donor im- 
purity, and the energy of the electron is strongly affected 
by quantum confinement within the NC. As an example, 
a single donor in the center of a NC has a delocalized elec- 
tron state when D < Gas [la 13 • This condition can be 
used as a somewhat conservative estimate for how small 
the diameter should be to produce electron states that 
are extended across the NC. 

In order to obtain the quantum energy spectrum in 



NCs, one can make the approximation that each NC is 
an infinite 3D square well. Such an approximation is valid 
because of the NCs' relatively large work function. The 
resulting energy spectrum can be described by defining 
the energy Eg (n) of the nth lowest electron, which gives 
for the first few energy levels 



Eqin) 



mD'^ 



0, n^O 

19.74, n = l,2 

40.38, 3<n<8 

66.43, 9 < n < 18 



(1) 



These first three nonzero energy levels can be labeled IS, 
IP, and ID, respectively. Higher electron shells have thus 
far not been examined by experiment, since they corre- 
spond to very large doping, and will not be discussed in 
this work. We focus primarily on the case where the spac- 
ing between quantum energy levels sa 20h^ /mD^ is larger 
than the characteristic scale of Coulomb energies, e^ /kD. 
The expression of Eq. ([1]) ignores the weak perturbation 
of quantum energy levels resulting from electron-electron 
interactions. This approximation is justified because of 
the large internal dielectric constant knCj a-s explained 
below. 

During the doping process, each NC i acquires some 
number Ni of positively-charged donors that it contains 
within its interior. These are assumed to be fixed, while 
the number of electrons ni within the NC can change 
due to electron tunneling between NCs. We assume that 
donors are added randomly to each NC by some high- 
temperature process, so that if the average number of 
donors per NC is v, then the probability that a given 
NC will have exactly TV donors is given by the Poisson 
distribution: 



,N 



^(^) - m'~ 



(2) 



This randomness in the number of donors is the only form 
of disorder that we include in our model. We show in 
Sec. IIVI that this disorder is sufficient to produce random 
charging of NCs, which leads to VRH. As mentioned in 
the introduction, the spontaneous charging of NCs is the 
result of the large gaps between quantum kinetic energy 
shells, which drive electrons away from NCs with many 
donors (emptying higher shells) and into NCs with few 
donors (filling lower shells), so that the number of elec- 
trons in a given NC is not generally equal to the number 
of donors. Additional disorder arising from fiuctuations 
in the NC size is not considered explicitly in this paper. 
The possible effect of such size fluctuations is discussed 
at the end of Sec. IIVI but we note here that fluctuation 
of NC size alone cannot produce spontaneous charging of 
NCs in the global ground state, which, as we show below, 
plays a crucial role for VRH. 

In addition to the quantum kinetic energy of the sys- 
tem, transport through the array is also greatly affected 
by long-ranged Coulomb interactions, which must be 
taken into account. In general, one could expect that 
calculating the total Coulomb energy of the system is 



a difficult problem, since the positions of negative elec- 
trons within each NC are described by their correspond- 
ing quantum wavefunctions and the positions of positive 
donors are random within the NCs volume. For our 
problem, however, a significant simplification is avail- 
able because the internal dielectric constant knc is much 
larger than both the external dielectric constant kj of the 
insulator in which the NCs are embedded and the overall 
effective dielectric constant k of the assembly. Specifi- 
cally, the large internal dielectric constant knc implies 
that any internal charge e is essentially completely com- 
pensated by the dielectric response, with the great major- 
ity of that charge, e(KNC ~ k)/knc7 becoming distribut- 
ing across the surface of the NC. In this way each NC 
can be thought of as metallic in terms of its Coulomb 
interactions. This allows us to write that the Coulomb 
self-energy of a NC with net charge q is given approxi- 
mately by q^ /kD, irrespective of how its constituent in- 
ternal charges are arranged. The interaction between 
two NCs i,j at a distance Vij can also be approximated 
as qtqjl KTij. These approximations are equivalent to the 
so-called constant interaction model, which is commonly 
used for individual quantum dots [17|. 

It should be noted that the effective dielectric constant 
K of the NC array is not simply equal to the dielectric 
constant ki of the insulating medium between NCs, but 
also includes the effect of polarization of NCs in response 
to an applied field. This polarization effectively decreases 
both the Coulomb self-energy of a single NC and the in- 
teraction between neighboring NCs. Generally speaking, 
the renormalization of the dielectric constant is not very 
strong, so that k is not very different from kj even when 
'^NC S> Kj. The canonical Maxwell- Garnett formula gives 
the approximate relation [18] 



(3) 



^ KNc + 2Ki -I- 2/(knc - ^i) 
' KNC + 2Ki - /(knc - Ki) 

where / — ■nD'^ /[Q{D')'^] is the volume fraction occupied 
by the NCs^or / < 0.4, this expression is accurate to 
within 8% |19|. As an example, for the case of a cubic 
lattice with D = 5 nm and D' = 6 nm (so that / = 0.3) 
and for KNc/'^^i = 5, one has k « 1.6Ki. 

Given this model, we can write down the Hamiltonian 
for our system as 



H 



E 



kD 



Y.^Q^^) 



k=0 



E 



e^{N,-n.,)iN,-n,) 



(4) 



Here, the first term describes the electrostatic self-energy 
of NC i, which has charge qi = e{Ni—ni), the second term 
describes the total quantum energy of the n^ electrons on 
NC i, and the last term indicates the Coulomb interaction 
between different NCs. 

The ground state for a particular system (a set of donor 
numbers {TVj}) is defined by the set of electron occupa- 
tion numbers {rii} that minimizes the Hamiltonian H. 



Given the ground state configuration, one can determine 

i f) 
the energy of the highest filled electron level, e[ , and the 

lowest empty electron level, e^\ at each NC i. Specifi- 
cally, 



-(/) 



Egin,) + 



e^[{N,^n,)^-{N,-n, + l)^] 



kD 



E 



e{Nj - Hj) 



... I^Tij 



(5) 



and 



ef^ =Eq{jm + 1) + 



^-[{N,-n,-lf-{N,-n,y 



kD 



E 



e{N,-n,) 



KTi 



(6) 



For the global ground state configuration, e/ < e^ for 
all i,j. As alluded to in the introduction, the definitions 
of £1 ' and ei'' describe only the lowest energy state of 
an electron or hole added to the site i. For this reason we 
refer to the density of states of these energy states e^ 
as the DOGS. 

The resistivity of the NC array is largely determined 
by the set of these ground state single-particle energies 
{e^ } and {e^ }• In the following section we show how 
these energy states can be used to calculate both the 
ground state electron DOGS g{e) and the resistivity p as 
a function of temperature and doping level. Note that in 
this problem every site is represented by two energies, in 
contrast to the canonical impurity band of lightly-doped 
semiconductors 20] , where every donor has only one rel- 
evant excitation energy. 

It is also important to note that in our model these 
donor electrons are assumed to be responsible for all con- 
duction. In other words, we assume that the tempera- 
ture T is low enough (and the doping level u is high 
enough) that donor electrons are much more abundant 
than electrons activated from the valence band. In prac- 
tical cases, this assumption is easily met: it requires 
only that the thermal energy fc^T be much smaller than 
the band gap energy Eg. More exactly, it requires that 
UbT < Eg/\n[K^cD^Eg/e'aBv''l\ 



III. COMPUTER MODELING 

In this section we describe our computational method 
for calculating the density of states and the resistivity 
at a given value oi v, T, and D. These calculations are 
based on a computer simulation of a finite, cubic array 
of L X L X L NCs, which proceeds as follows. First, we 
specify the doping level v. The simulation then assigns 
the donor number Ni for each NC i according to Eq. ([2]) . 
The initial values of the electron numbers {rii} are then 
assigned randomly in such a way that the system is over- 



then searches for the ground state by looping over all NC 
pairs (ij) and attempting to move one electron from i to 
j. If the move lowers the Hamiltonian H, then it is ac- 
cepted, otherwise it is rejected. Equivalcntly, one can say 
that for each pair i,j we check that two ES ground state 
criteria are satisfied: 



and 
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KTi 
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KTi 
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(7) 



(8) 



all electro- neutral, i.e., J2i 



J2i^j- The simulation 



If either one of these criteria is violated, then an electron 
is transferred. This process continues until all sites i,j 
satisfy Eqs. (0 and ©. 

It should be noted that this procedure does not in gen- 
eral find the exact ground state, but only a "pseudo- 
ground state" that is stable with respect to single- 
electron transfers. In principle, the system energy can 
be lowered further by some multi-electron transfers. The 
effect of these higher-order relaxation processes on the 
properties of the pseudo-ground state has been exam- 
ined for similar models [211, ^^ , and they are generally 
beyond our intended accuracy in this paper, so we do not 
consider them here. 

Once the pseudo-ground state occupation numbers 
{rii} have been found, one can define the single-particle 

energies e/ and s^ for each NC i using Eqs. ([5|) and 
([5]). These energies are tabulated and then histogrammed 
in order to calculate the single-particle DOGS g{e). In 
the results presented below we define electron energies 
e relative to the Fermi level /i, which is calculated for 
each realization of the simulation as /i = [minJEj } — 
max{ei }]/2. In this way e < corresponds to filled 
electron states e^-^^ while e > corresponds to empty 
states e^'^\ (See, for example. Fig. 2] below.) 

Once the pseudo-ground state energies {e/ } and 

{el*^ } are determined, we calculate the resistivity of the 
system by mapping the simulated NC array to an ef- 
fective resistor network. The equivalent resistance Rij 
between NCs i and j can be determined by writing down 
the time-averaged rate of electron transfer between sites 
i and j in the presence of an electric field and expand- 
ing in the limit of small field, as in the canonical Miller- 
Abrahams resistor network [20, [23| . In calculating Rij we 
consider only electron transfer among the highest filled 
states, e*^-''), and the lowest empty states, e''^-', which is 
appropriate when the temperature is small enough that 
T < e^ /kD, so that thermal excitation of multi-electron 
transitions is exponentially unlikely. 

Since each NC has two energy levels that can partici- 
pate in conduction, e^-'') and e^''\ one can say that there 
are four parallel conduction processes that contribute to 
the resistivity between two NCs i and j: one for each 
combination of the initial energy level at site i (either e^ 



or e^^-^) and the final energy level at site j (either ey or 

e]j^ ). Each of these four processes has a corresponding 

effective resistance Rl" , where a,/3 = (/), (e). These 
four resistances can be said to be connected in parallel 
between NCs i and j, and their value can be written 
compactly as 



r'^^^ = Ro exp 



2r^ 



ij 

knT 



(9) 



where Rq is a prefactor that has only a relatively weak 
power-law dependence on temperature. The first term 
in the exponential of Eq. © describes the exponential 
suppression of the tunneling rate with distance r, as 
explained in the introduction, and the second term de- 
scribes thermal activation by exponentially-rare phonons 



of energy e 



(",^) 



Since we are interested only in iden- 



tifying the exponential component of the dependence of 
resistivity on temperature, we take Rq to be a constant. 
The energy e"' in Eq. © is defined as follows j20| : 



-(/3) 



-(",« 



-Wi 
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ef )£(") < 
efe^^ > 



(10) 



The net resistance Rij between NCs i and j is the par- 



allel sum of the four resistances R 



(a/3) 



Since the expo- 



nential factor in Eq. ([9]) provides a sharp differentiation 
between these four parallel resistances, at relatively low 
temperatures and to within the accuracy of our calcula- 
tions we can equate Rij with the minimum of the four 
parallel resistances. That is, 



Ri 



^in{i?r^} 



(11) 



After calculating all resistances Rij for a given simulated 
array, we find the dimensionless resistivity of the network 
p/poi where po = RqD' , using a percolation approach 
[20|. Specifically, we find the minimum value Re such 
that if all resistances Rij with Rij < Re are left intact 
while others are eliminated (replaced by Rij = oo), then 
there exists a pathway connecting the left and right faces 
of the simulation volume (the "infinite" percolation clus- 
ter). The resistivity p/po is approximated as Rc/Rq- 

In our analysis below we make use of the following di- 
mensionless units, which reduce the number of free vari- 
ables in the problem. We introduce the dimensionless 
distance between the centers of NCs i and j. 






the dimensionless temperature 

2D^KkBT 



rp^ 



e^i 



(12) 



(13) 



the dimensionless electron energy 

e 



the dimensionless electron DOGS 

5*(e*) = g(£*), 



(14) 



(15) 



and the dimensionless resistivity 



Inp* 



2D 



ln(p/po)- 



(16) 



In these units, Eq. (J9l) can be written more simply as 



In pI 



JT* 



(17) 



and the problem loses any explicit dependence on the 
diameter or the localization length. It is also convenient 
to discuss the energy gap between the IS and IP shells 
in terms of the dimensionless parameter 



EQ{i)^EQ{2) 
e^/nD 



20.64 



me^D 



20.64 



kqb 



(18) 

We use our simulation to examine the resistivity at vari- 
ous values of ly, T*, and A. 

Results below correspond to a simulated system of size 
L — 25 with open boundaries, averaged over 100 realiza- 
tions. Simulations at smaller system size, 15 < L < 25, 
do not produce noticeably different results for either the 
DOGS or the resistivity, which allows us to avoid having 
to extrapolate our results to infinite system size. 



IV. RESULTS AND DISCUSSION 

Our goal is to determine which conditions produce 
VRH in the NC array. To this end we calculated the 
resistivity p and the electron density of ground states 
g{e) for a range of values of the doping level i/, the tem- 
perature T* , and the quantum energy scale A. (Varying 
A is equivalent to considering different values of the NC 
diameter.) Before proceeding to present general results, 
however, we first illustrate the most important features 
of the problem by discussing the hypothetical case where 
all NCs have the same number of donors, so that there is 
absolutely no disorder in the system. Say, for example, 
that 1/ = 5 and that Ni = 5 for all i. In this situation, 
the ground state arrangement of the system is for elec- 
trons to uniformly neutralize all donors: rn = Ni — 5. 
The result, by Eqs. ([5]) and |6]), is that every NC has 
the same two energy levels, e'^^^ — Eq{5) — c^/kD and 
e(«=) — Eq{5) + e^/nD, and the system's Fermi level 
/i — Eq{5). Equivalently, one can say that the single- 
particle DOGS for this hypothetical system corresponds 
to two (5-function peaks at e = ±e^ / kD. 

As explained in the introduction, conduction in this 
uniformly neutral system requires the excitation of a pos- 
itive/negative NC pair. Specifically, such an excitation 



produces one positive NC containing 4 electrons and one 
negative NC containing 6, and as such it has an excitation 
energy equal to the sum of the two Coulomb self-energies. 
Equivalently, one can say that conduction requires the 
production of a hole in the filled (5-function DOGS peak 
at e = ~e^ /kD and an electron in the empty DOGS 
peak at e = e^ / nD, and so the conduction has an acti- 



vation energy ea 



' /kD. Thus, this hypothet- 



ical system without disorder has activated conduction: 
p^ PQeyip[eA/kBT]. 

On the other hand, once the randomness in donor num- 
ber is taken into account, one can no longer say in gen- 
eral that the ground state arrangement of electrons is 
uniformly neutral, rii = N^. Indeed, when N^ can take a 
wide range of values, then those NCs with very large N 
may become ionized so that their electrons can occupy 
lower-energy shells on other NCs with small N. In this 
way, the presence of a discrete quantum energy spectrum 
instigates the production of positively- and negatively- 
charged NCs. It is this spontaneous charging that allows 
for VRH, as we will show below. 

Still, it is straightforward to see that the system re- 
mains nearly uniformly electroneutral in the ground state 
under either of two conditions: (i) very small quantum 
energy gap, A ^ 1, or (ii) very small doping level, j^ ^ 1. 
In the former case, the difference between quantum en- 
ergy levels becomes negligibly small compared to the en- 
ergy required to produce charging of NCs. Thus, the 
NCs remain neutral and the conduction is activated, as 
explained above. In the limit of very small doping, ly <^ 1, 
the system also remains nearly uniformly neutral due to 
an extreme scarcity of donors with Ni > 2. Indeed, by 
Eq. ^, at small v the fraction of donors with Ni > 2 
is ~ ly'^ /6. Thus, neutrality of the system can be main- 
tained without requiring any significant number of elec- 
trons to occupy the IP shell, and there is essentially no 
charging of NCs. Therefore in the limit of very small i/ 
the conduction is also activated. 

In situations where either A or j/ is not small, one 
can expect spontaneous charging of NCs in the ground 
state, and it is not trivial to predict the DOGS or the 
temperature dependence of the resistivity. We explore 
these situations using our simulation method, outlined in 
Sec. mil Before proceeding to present results for a wide 
range of ly and A, we first focus on the illustrative cases 
oi v — 5 and v = 2, taking for the quantum energy gap 
A = 5. 

At 1/ = 5, the Fermi level resides in the middle of the 
IP shell. Thus, since the gap between quantum energy 
levels is relatively large, in the ground state essentially 
all NCs satisfy 2 < n^ < 8. By Eq. ([2]), however, roughly 
11% of NCs have a donor number satisfying A^i < 2 or 
Ni > 8. Such NCs become charged in the ground state, 
driven by the large gaps in the quantum energy spectrum 
that induce electrons to leave the ID shell and to fill the 
IS shell. Thus, the ground state configuration of the 
system consists of randomly-distributed fixed charges, 
which correspond to those NCs with Ni < 2 (which be- 



come negatively-charged) or Ni > 8 (positively-charged) , 
and mobile electrons and holes in the partially-filled IP 
shell. The mobile electrons and holes arrange themselves 
in such a way that the ES criteria of Eqs. (O and (|S]) are 
satisfied. It is these criteria that giv e rise to the vanishing 

iiiia. 



DOGS near the Fermi level 



This process of charging of NCs is illustrated schemat- 
ically in Fig. [21 which shows the energy levels of isolated 
NCs with donor numbers < iV < 10. In the neutral 
state, a NC with N donors has N filled electron energy 
levels (Fig. [3^) . When the system contains a mixture of 
NCs with different N, however, electrons abandon high 
energy levels in NCs with large N and fill empty states 
in NCs with small N. This process is shown for the case 
ly = 5 in Fig. [3)d. The resulting charged NCs produce a 
random Coulomb potential throughout the system that 
smears the single electron energy levels and produces a 
finite density of states near the Fermi level. 

The DOGS for j/ = 5 and A = 5, as calculated by 
our numerical simulation, is plotted in Fig. 3^. One can 
see the quadratic Coulomb gap near the Fermi level, as 
proscribed by the ES theory. As compared to the con- 
ventional Coulomb gap problem in lightly-doped semi- 
conductors ^20,] , this Coulomb gap is remarkably well pre- 
served, with the DOGS remaining quadratic until e* « 1. 
This strong Coulomb gap suggests that the resistivity 
should follow the ES law for all temperatures T* ^ 1. 
Specifically, at these small temperatures the resistivity is 
described by 



PiT) == Po exp 



T, 



ES 



T 



1/2- 



where 



'^ES 






(19) 



(20) 



and C is a numerical coefficient of order unity. 

This behavior can indeed be seen in Fig. [Jb, where 
Inp* is plotted as a function of (T*)~^/^. The linear re- 
lationship at large (r*)~i/2 suggests that, as expected, 
the resistance follows the ES law at small temperatures. 
We find that the numerical coefficient C ~ 8.1, as com- 
pared to the typical value C ~ 2.8 in lightly-doped bulk 
semiconductors |20i]- At larger temperatures T* > 1 [or 
(y*)-i/2 <; 1]^ the resistivity saturates at Inp* = 1. At 
such large temperatures the factor e*JT* in Eq. pT]) typ- 
ically becomes smaller than unity, which indicates that 
electrons tunnel relatively easily between nearest neigh- 
bors, and VRH is abandoned in favor of nearest- neighbor 
hopping. At these large temperatures the resistivity can 
be expected to have only a relatively weak power-law de- 
pendence on temperature, which is beyond the accuracy 
of our numerical calculations. 

In addition to the parabolic Coulomb gap near the 
Fermi level, another salient feature of the DOGS in Fig. 
S^ is that it has strong maxima at e* — ±1 and col- 
lapses nearly to zero at e* = ±2, as if there were addi- 
tional Coulomb gaps that constrain the density of states 



a) 



v = 5 



N = 123456789 10 



2+A 



£it-T- 



b) 




P(N) , % : 0.7 3 



18 15 10 



FIG. 3. (Color online) Schematic depiction of the charging 
process in a system with NCs with varying donor number 
A^. a) The single-electron energy levels (horizontal line seg- 
ments) are shown for isolated NCs. The Coulomb self-energy 
of charged NCs produces a spectrum where different charge 
states have a separation 2e^ /nD. The quantum confinement 
energy provides a gap between subsequent shells, e.g. IS and 
IP states or IP and ID states. In the neutral state, a NC 
with A'^ donors has A'' filled energy levels (indicated by filled 
blue dots), els indicates the quantum kinetic energy of the IS 
shell, £is ~ ^q(1)/(6^/'^^)- t)) A depiction of the charging 
process at i^ = 5. Electrons in the ID shell of NCs with N > 8 
abandon these NCs and instead fill empty energy levels in the 
IS shell of NCs with A^ < 2. In this way NCs with A^ > 8 
become positively charged and NCs with N < 2 become neg- 
atively charged. The resulting Fermi level /x is shown by the 
dashed line. For NCs with N = 5, it resides in the center 
of the IP shell. The relative abundance of different donor 
numbers at i^ = 5 is shown at the bottom of the figure as a 
percentage. 



around e* ~ ±2. These "reflected Coulomb gaps" are 
in fact the product of an approximate symmetry in the 
system, which can be seen by examining Eqs. ([5]) and ([6|). 
At z^ = 5, the great majority of NCs have 2 < n^ < 8. 
For such NCs, Eg^ni) — Eglrii + 1); both the high- 
est filled and lowest empty electron states are in the IP 
shell. In this case, one can subtract Eqs. ([5]) and ^ to 



show that ej = e. 



K/) 



2. Thus, the great majority of 
NCs contribute to the density of states two energy lev- 
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FIG. 4. (Color online) Density of ground states and resistivity 
at i^ = 5 and A = 5, as measured by computer simulation, a) 
Density of states as a function of electron energy. Filled elec- 
tron states are shaded. Dashed red lines show, schematically, 
the quadratic Coulomb gap near the Fermi level, e* — 0, and 
the "reflected Coulomb gaps" at e* = ±2. Note that that the 
total shaded and unshaded areas under the g*{e*) curve are 
both normalized to unity, since each NC has one electron and 
one hole excitation. The inset shows the DOGS over a wider 
energy range, with small, distant peaks indicating rare NCs 
whose highest filled electron state is in the IS shell or whose 
first empty state is in the ID shell, b) The dimensionless 
logarithm of the resistance, Inp*, as a function of (T*)"^", 
which illustrates the existence of ES resistivity at small tem- 
perature. 



els - one filled, one empty - separated by 2e'^ / kD. This 
creates an approximate discrete translational symmetry 
in the density of states, so that g*{e*) « g*{e* - 2) for 
< £* < 2. As a consequence, the Coulomb gap at the 
Fermi level implies the existence of reflected Coulomb 
gaps at e* — ±2. In other words, one can say that be- 
cause of the discrete charging spectrum of NCs the con- 
ventional quadratic bound on the DOGS near the Fermi 
level also produces (approximate) quadratic bounds on 
the DOGS near e* = ±2. The contribution of rare NCs 
with rii — 2 OY rii — % to the DOGS can be seen in the 
small peaks at e* = —6 and e* — 7, as shown in the inset 
of Fig. Hi. 

The presence of reflected Coulomb gaps is not unique 
to the doping level v = ^. Indeed, for all v that are 
sufficiently removed from the quantum energy gaps at 



_*(e) 



e*^^^ + 2 is valid 



1/ = 2, 1/ = 8, etc., the relation e 
for most NCs in the system and the resulting DOGS is 
essentially identical to that of Fig. \^. Consequently, 
the resistivity plot shown in Fig. HJd accurately describes 
the resistivity at most values of z/ > 1. The reflected 
Coulomb gaps in Fig. |4^ appear even more dramatically 
for large NCs with external impurity charges, as will be 
shown in Sec. |Vl 

On the other hand, one could expect qualitatively 
different behavior at z^ = 2, where there are precisely 
enough electrons to fill the IS shell of every NC, and the 
Fermi level sits in between the IS and IP shells. In this 
case there is no "discrete translational symmetry" in the 
density of states, since the empty and filled energy levels 
for most NCs, e] and e] , sit on opposite sides of the 
quantum energy gap, as shown schematically in Fig. [5] 
This produces a DOGS that is qualitatively different from 
what is shown in Fig. |4^. One could therefore expect 
that the dependence of the resistivity on temperature is 
also qualitatively different. Such thinking is supported 
by a recent experiment on electrochemically gated NCs 
[3|, which reported that when v is very close to 2 there 
appears an appreciable temperature window over which 
the resistivity follows the Mott law. Given these differ- 
ences, it is worth giving some special consideration to the 
case 1/^2. 



IS and IP electron states contribute to the DOGS near 
the Fermi level, the density of states has a characteristic 
width similar to that of the gap A. 
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FIG. 5. (Color online) Schematic depiction of the filled and 
empty energy levels at // = 2. Energy levels are shown for 
NCs in the absence of any Coulomb potential, similar to Fig. 
|3l At z^ = 2, some electrons leave the IP shell of NCs with 
A'' > 2 and fill empty states in the IS shell of NCs with A'' < 2. 
The resulting Fermi level /x is aligned with the first(second) 
energy level of the IP shell in NCs with A'' = 4(5), which is 
partially filled. 

The DOGS for ;/ = 2 is shown in Fig. Hi. Unlike at 
i^ = 5, where the DOGS collapses at e* = ±2, the DOGS 
at i^ = 2 is much broader, with a width A + 2. This 
broad DOGS can be seen as a consequence of the large 
gap between IS and IP energy shells, which implies that 
the energy of electron or hole excitations, e] ' and e^', 
can take a wide range of values, depending on the donor 
number iV^. Alternatively, one can say that since both 



FIG. 6. (Color online) Density of states and resistivity at 
V = 2 and A = 5, as measured by computer simulation, a) 
DOGS as a function of electron energy. Filled electron states 
are shaded. The dashed red curve is the same parabolic curve 
shown in Fig. |4^. The inset shows the DOGS very close to the 
Fermi level, b) The dimensionless logarithm of the resistance, 
Inp*, as a function of (T*)"^", which shows ES resistivity at 

As at 1/ = 5, the DOGS vanishes at the Fermi level 
(see the inset of Fig. |6K), but in this case it can only 
be described as parabolic over the fairly narrow range of 
energies |e*| < 0.2. In the intermediate range of energies 
0.2 < |e*| < 1, the DOGS grows roughly linearly with en- 
ergy. At larger energies 1 < e* < A the DOGS becomes 
roughly constant. 

In spite of this relatively comphcated DOGS, Fig. |61d 
shows that the resistivity is in excellent agreement with 
the ES law, with a coefficient C w 5.7 [see Eq. pi)) ], 
at all but very large temperatures. This is somewhat 
surprising, since it suggests that the system exhibits ES 
resistivity even when the temperature is large enough 
that the band of energies over which VRH occurs is much 
larger than the width of the parabolic Coulomb gap. This 
behavior would be impossible if states were randomly 
distributed in space. Our observation of ES resistivity 
suggests that at z/ = 2 spatial correlations emerge which 
somehow preserve ES resistivity even in the absence of a 
parabolic DOGS. 

To illustrate how this might be possible, let us first 
recall that in a disordered two-dimensional (2D) system, 
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the DOGS is linear in energy near the Fermi level rather 
than parabolic, but the ES law of VRH is still obeyed 
[12| . One can now imagine a 3D system in which sites 
with energies close to the Fermi level are arranged in 
a 2D fractal subspace embedded in the system volume. 
In such a system, one would still have a linear DOGS 
near the Fermi level accompanied by ES resistivity, even 
though the system as a whole is three-dimensional. Using 
this reasoning, one can speculate that the results shown 
in Fig. ini are indicative of such a fractal arrangement of 
sites near the Fermi level, driven in some way by the long- 
ranged Coulomb potential. More broadly, these results 
hint at the idea that in a disordered system of localized 
states dominated by Coulomb interactions, one should be 
able to derive the ES law without explicit reference to the 
DOGS or the system's dimensionality. Such an argument 
was in fact first put forward by Larkin and Khmelnitskii 
[2J|. Our system at :/ = 2 may be a good application 
of this argument. It remains unclear, however, in which 
situations this argument is applicable a priori. This gen- 
eral question and its application to the case v = 2 will 
be the subject of a future publication. 

By conventional thinking, the relatively constant 
DOGS at \e*\ > 1 would seem to suggest a regime of 
temperature in which the resistivity follows the Mott 
law, which describes VRH in the presence of a constant 
DOGS. However, unlike the experiments of Ref. [3], we 
see no noticeable region of Mott VRH. The Mott resis- 
tivity observed in Ref. Q at i^ = 2 is likely the result of 
some additional disorder that is outside the model con- 
sidered in this section, and is discussed further in Sec. 

m 

Having considered the specific cases of t^ = 5 and 
V = 2, we now turn our attention to a general description 
of VRH at different values of i' and A. In order to iden- 
tify more precisely which conditions produce VRH, we 
used our simulation to measure the resistivity as a func- 
tion of T*, v, and A over the range 0.01 < T* < 10, 
0.2 < t^ < 2, and 0.5 < A < 5. For each case we 
measured the exponent 7 of the temperature dependence 
of resistivity by calculating the "reduced activation en- 
ergy" w{T*) = -d{lnp*)/d{\nT*) oc T-t ^. The ex- 
ponent 7 was identified by making a power law best fit 
to w{T*). Those values of T*, v, and A that produce 
7 — 0.5 ±0.1 were identified with ES resistivity; domains 
where 7 > 0.6 were identified with activated resistivity. 
As discussed above, no significant regimes were identified 
that showed Mott behavior. We use this data to con- 
struct an approximate phase diagram in the space of T* , 
V, and A that identifies which behavior can be expected. 

Our result is plotted in Fig. [2] for T* ^ 1. Gener- 
ally speaking, the results indicate that for v > 0.6 and 
A > 0.5 one can expect ES resistivity, while for other 
conditions the resistivity is activated. These conditions 
are equivalent to the conditions (i) and (ii) that were 
announced in the introduction. Dashed horizontal lines 
indicate, as an example, the values of A corresponding 
to CdSe NCs with D = 6.2 nm, as in Ref. |7], and Si NCs 



with D = 5 nm, as in Ref. [l| . Both of these dashed lines 
assume that knc/^^i — 5- At temperatures T* > 1 VRH 
is gradually replaced by nearest-neighbor hopping. The 
condition T* < 1 is equivalent to the condition (iii) from 
the introduction. 

As mentioned above, the model considered in this sec- 
tion does not account explicitly for any sources of disor- 
der other than fluctuations in donor number. For exam- 
ple, in real NC arrays the diameter D varies from one NC 
to another, which introduces variations in the quantum 
spectrum between NCs [see Eq. ([T])]. Nonetheless, the 
presence of these size fluctuations in addition to fluctua- 
tions in donor number does not destroy ES VRH, since 
the Coulomb gap near the Fermi level is a universal result 
of the ES stability criteria [Eqs. ([7]) and ^] and is inde- 
pendent of the source of disorder in the system. Whether 
size fluctuations or other sources of disorder enhance the 
role of VRH or significantly affect the magnitude of the 
resistivity remains yet to be studied. Generally speaking, 
however, one can expect that the phase diagram of Fig. 
[5] is accurate whenever the typical magnitude of size fluc- 
tuations 6D satisfies {SD)/D ^ 1/A. Wc further expect 
that even larger size fluctuations do not greatly affect 
VRH in regimes where the ES law applies, since in such 
cases the DOGS is already saturated by the disorder in 
donor number. In regimes where the resistivity is acti- 
vated, the presence of a large additional disorder should 
generally promote the existence of VRH, which decreases 
the resistivity at small T*. 



V. VRH IN ARRAYS OF LARGE NCS WITH 
EXTERNAL IMPURITY CHARGES 

In the previous sections we showed that the sponta- 
neous charging that leads to VRH is driven by the rela- 
tively large gaps between degenerate shells of the electron 
quantum energy spectrum. In large NCs, the gap A be- 
comes small and this charging disappears, which leads 
to activated resistivity (see Fig. [5]). A similar effect can 
be expected when the NC shape is not symmetric. In 
this case, the electron energy levels are not degenerate, 
so that the bunch of energy levels corresponding to a 
particular shell in a spherical NC is dispersed, and as 
a consequence the gaps between subsequent energy lev- 
els are reduced. Thus, large or highly asymmetric NCs 
tend to remain neutral in the ground state and exhibit 
activated transport even at i^ > 1. 

In this section, however, we show that if donor impu- 
rities are located outside of NCs, ES VRH can still be 
observed. The presence of ES VRH in large NCs can be 
understood using an argument that was first put forward 
by Ref. 10 in the context of granular metallic films. The 
argument from Ref. |l0| is briefly repeated here. 

Consider an array of spherical semiconductor NCs with 
large internal dielectric constant knCj each of which is 
coated with a thin layer of width w of insulator, such as 
the semiconductor's own oxide or the ligands shown in 
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Fig. [TJ Suppose further that donor impurity charges +e 
are embedded in this insulator, as shown schematically 
in Fig. [71 with some overall concentration iVimp- If some 
particular donor resides at a location within the insulator 
shell that is well-separated from the points of contact be- 
tween neighboring NCs, then this donor simply donates 
its electron to the NC on which it resides. The resulting 
positive impurity charge induces a negative image charge 
on the NC surface because of the dielectric response, and 
together the donor and its image charge make a com- 
pact neutral pair. In this way, donors that are not near 
the point of contact between two NCs produce a negligi- 
ble Coulomb potential that plays no role in charging the 
system. 

On the other hand, when a donor is located close to 
the point of contact between two NCs, labeled A and 
B, it induces negative image charges —qA and — g^ in 
the surfaces of NCs A and B, respectively. In order to 
maintain overall neutrality of the NCs, an equal and op- 
posite image charge appears at the center of each NC: 
+qA and +qB- (These "image charges at the center" 
represent a uniform electronic charge at the NC surface.) 
The values of qA and qB are such that together the im- 
age charges —qA and —qs neutralize the donor charge: 
qA+qs = e. Their respective magnitudes are determined 
by the distance between the impurity and each NC sur- 
face. For example, if the impurity sits exactly along the 
line connecting the centers of NCs A and B and if the 
gap 2w between NCs satisfies w -^ D, then the NCs 
can be approximated as infinite planar metallic surfaces 
and qAXB — qs^A, where xa and xb are the distances 
between the impurity and the surface of NCs A and B, 
respectively. More generally, when w <^ D, the image 
charges —qA and —qB sit very close to the donor im- 
purity and essentially neutralize it, so that the screened 
impurity does not directly contribute to any Coulomb 
potential at length scales of D or larger. Instead, the net 
effect of the image charges is to "fractionalize" the donor 
impurity charge, such that +qA is relayed to the center of 
NC A and +qB is relayed to the center of B. This process 
is depicted in Fig. [T] 

In this way, each NC gets a number of random, positive 
fractional donor charges created by those donors located 
near the contact points between NCs. We denote the 
sum of all fractional charges at NC i by Qi . The propor- 
tion of all donor charges that sit at these contact points is 
~ w/D, so that of the total number r^ Nj^^pwD^ of donor 
impurities in the insulator shell covering a given NC, only 
~ NimpW^D of these become fractionalized. When the 
average total number of fractionalized charges per NC 
NiaipW^D :$> 1, the central limit theorem guarantees that 
the distribution of the random variable Qi/e can be ap- 
proximated as a Gaussian with mean ~ Ni^ipW^D and 
root mean square fluctuation r^ {Ni^pw'^Dy/^. 

Donor electrons respond to the potential created by 
fractional charges by arranging themselves on NCs in in- 
teger number and in such a way that the total electro- 
static energy of the system is minimized. In other words. 
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FIG. 7. (Color online) A schematic depiction of the fraction- 
alization of the charge of a donor impurity between large NCs 
(see also Fig. 11 of Ref. 10|). Semiconductor NCs (tan/gray 
circles) have a thin coating of insulator (yellow/light gray), 
with embedded donor impurities (small black circles). Each 
positive donor induces negative image charges (small white 
circles) that neutralize it, while equal and opposite positive 
images are conveyed to the center of the NC. Those donors 
that are located close to the point of contact between two 
NCs create non-integer image charges in the two surrounding 
NC surfaces. In this way NCs are given a "fractional donor 
charge" Qt. 



in the ground state the set of electron occupation num- 
bers {rii} is that which minimizes the Hamiltonian 



^-E 



{Qi - euiY 
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Unlike in the Hamiltonian of Eq. ^ of Sec. HIl here the 
quantum energy gaps are negligibly small and disorder is 
provided by the fractional charges {Qi}- 

Eq. (PT|) implies that the corresponding electron energy 
levels at NC i are given by 

(/) _ [Qi - eriiY - {Qi - erij + e)^ ^^ Q-j - erij 



and 



(e) _ iQt_ 



kD 



ef - (Q. 



... Hfij 



(22) 



g, 



kD 



i/j 



(23) 

for the highest filled and lowest empty states, respec- 
tively. By subtracting Eqs. (1^^ and (1^^ it can be seen 
that e^"'^ = e^f^ +26^ / kD for all NCs. This has important 
implications for the DOGS, as will be shown below. 

As explained above, the values of the fractional charges 
Qi can be assigned using a Gaussian distribution with 
mean Q ^ eNi^pW^D and standard deviation uq = 
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{eQY/"^. One can notice, however, that in the Hamil- 
tonian of Eq. (|21l) the variables Qi and rii appear only 
in the combination Qi — eui, which by electroneutrality 
of the system must satisfy {Qi — erii) — 0. Thus, when 
calculating the DOGS and resistivity, one can adopt a 
somewhat simpler model where Qi is chosen from a dis- 
tribution with mean zero and rii is allowed to take any 
integer value (positive or negative). 

In fact, an even further simplification of the model 
is available when the standard deviation ctq/b ^ 1. 
Namely, when Qi can take such a wide range of values, 
one can approximately replace the broad distribution for 
Qi/e with a uniform distribution Qi/e G [—1/2, +1/2]. 
The validity of this approximation can be understood by 
considering that each NC minimizes its Coulomb self- 
energy by minimizing the magnitude of its net charge, 
\Qi — eni\. Since Ui can take any integer value, it is gener- 
ally true that in the ground state — e/2 < Qi — erii < e/2. 
This random spatial arrangement of net charges produces 
a fluctuating Coulomb potential that leads to ES VRH 
[26| . All results below correspond to this choice of a uni- 
form distribution for Qi. 

It is worth noting that while so far we have focused 
on positive donor impurities, acceptors embedded in the 
insulating layer can play the same role. Random charging 
of NCs can also result from the simultaneous presence of 
both donors and acceptors. For the sake of argument, 
however, we focus our discussion around positive donor 
charges. 

Results for the DOGS and resistivity are given in Fig. 
[SI as calculated using the simulation method described 
in Sec. IIIII and the definition of single-particle energies 
given in Eqs. ^^ and ((231). One can note that the DOGS 
shown in Fig. |S^ vanishes at the Fermi level e* — 0, as 
required by the ES criteria [Eqs. ^ and (|S])], as well as 
at e* — ±2. Because of the lack of a quantum energy 
term in the Hamiltonian, the DOGS also has a perfect 
discrete translational symmetry: g*{e*) = g*{s — 2) for 
e* > 0. The strong "reflected Coulomb gaps" in Fig. [5] 
are a result of the relation e* = e* +2. 

The resistivity for this system is plotted in Fig.[B]3 as a 
function of (T*)~^/^. As expected, the resistivity follows 
the ES law at T* ^ 1 due to the strongly-preserved 
quadratic Coulomb gap near the Fermi level. The dashed 
line in Fig.[5]D corresponds to the ES law with a coefficient 
C«9.6 [seeEq. ^]. 

We have also verified that our results for the DOGS 
and resistivity are practically identical if Qi/e is chosen 
not from a uniform distribution [—1/2, -1-1/2] but from a 
Gaussian distribution with three-times larger variance. 

In the opposite limit, where fractionalized donor 
charges are very rare, Nij^pW^D ^ 1, each NC remains 
essentially neutral, and there is no random Coulomb 
potential. This uniformity leads to activated nearest- 
neighbor hopping, since without disorder long-range elec- 
tron hops cannot reduce the energy required for the hop 

Finally, it can be noted that in our discussion above 
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FIG. 8. (Color online) Electron DOGS and resistivity for 
an array of large semiconductor NCs with fractional donor 
charges corresponding to Qi/e £ [-1/2,-1-1/2]. a) The 
DOGS, which vanishes at e* = and at e* — ±2 because 
of the ES stability criteria and the relation between filled and 
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which shows ES behavior at temperature 



we have ignored the possible presence of deep electronic 
states at the NC surface. Such trap states can play the 
role of compensating impurity centers, which remove a 
percentage of electrons from the conduction band. In 
this case, only the uncompensated donor electrons con- 
tribute to conduction, and the value of v is effectively 
renormalized downward. Repulsion between uncompen- 
sated donor electrons and electrons in trap states may 
also produce a small shift in electron energies, and is 
outside the scope of our treatment here. 



VI. GATING OF A NC ARRAY BY AN IONIC 
LIQUID 

In Sees. |T]-|V]we discussed systems of NCs doped by 
random impurities, and we explored the dependence of 
the resistivity on the doping level. In such systems, the 
doping level is established during the fabrication of NCs. 
In many cases, however, it is desirable to have a doping 
level that can be continuously tuned, so that the resis- 
tivity of a single device can be set to a wide range of 
values. For this purpose, electrochemically gated arrays 
of semiconductor NCs are actively being studied 0, 0] ■ 

In such systems, conduction electrons are introduced 
into the system via a voltage source, which drives elec- 
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trons from a top gate to a bottom gate that is in electrical 
contact with the NC array. Generally, in between the top 
gate and the NC array is a room temperature ionic liq- 
uid that provides large capacitance and therefore allows 
for a high density of electrons to be introduced to the 
NC array at a relatively small voltage [28|. The cations 
from this ionic liquid intercalate into the spaces between 
NCs, penetrating deep into the array through the perco- 
lating network of pores between NCs, and thus provide 
a neutralizing charge for the conduction electrons. A 
schematic picture of this system is given in Fig. [5] 




FIG. 9. (Color online) A schematic picture of an array of 
semiconductor NCs (large circles) gated by an ionic liquid. 
Cations (small circles with -|-'s) are driven by a voltage source 
to intercalate between NCs. Because of the large NC dielectric 
constant knc, the net effect of positive ions is to provide a 
fractional donor charge Qt at a given NC i, similar to what is 
shown in Fig. [7] Neutralizing electrons occupy NCs in order 
to neutralize ionic charges. Ligands separating NCs are shown 
as curvy lines. 

The large internal dielectric constant of NCs and the 
relatively small diameter of cations suggests the presence 
of strong image charge forces that bind cations electro- 
statically to their image charges in the NC surface. In 
this way, one can expect that cationic charges are located 
primarily on the surface of each NC. If one assumes that 
the position of cations on the NC surfaces is random, 
then one again arrives at a model of fractionalized cation 
image charges, similar to what is suggested in Sec. El 

For this model one can use a Hamiltonian that includes 
both a prominent quantum kinetic energy spectrum, as 
in Sec. |lTl and a fluctuating, fractionalized donor charge, 
as in Sec. El 



^-E 



[Qi - eriiY 



kD 



E^q(^) 



E 



fc=0 

(Qt - eni){Qj - en.,-) 



(24) 



Here, the fractional charge Qi/e can be chosen uniformly 
from the interval \v — 1/2, v + 1/2]. 

Using our computer simulation method, we have briefly 
investigated the DOGS and resistivity of the system de- 
scribed by this Hamiltonian at various values of i^ > 1. 
We find that ES VRH appears at low temperature for all 
values oi V > \. In fact, when |z^ — 2| > 1 and ji' — 8| > 1, 
the DOGS is exactly the same as in Fig. [S^, and the 
resistivity is also identical. 

We note that the model defined by Eq. ([24]) . where the 
fractional donor charge is completely random, is unlikely 
to be accurate when v is at the boundary between two 
quantum energy shells. At z/ = 2, for example, random 
fractional charges lead to a fluctuating Coulomb potential 
with characteristic amplitude much larger than ksT /e at 
room temperature. However, such a large Coulomb po- 
tential induces cations, which are mobile during the gat- 
ing process, to rearrange in order to screen the potential. 
In this way the cation positions become correlated and 
the typical amplitude of the Coulomb potential is reduced 
to fcsT/e, which is not large enough to produce charging 
of NCs. As a result, the typical amplitude of fluctuations 
in Qi is likely much smaller than e, so that one should not 
expect a finite DOGS near the Fermi level. Rather, in 
the absence of any other disorder, the resistivity should 
be large and activated. 

Experiments with ionic liquid gating confirm that, as 
expected, the resistivity is much larger at t^ = 2 than at 
other filling factors [7|. However, the resistivity is gener- 
ally shown to correspond to VRH rather than activated 
behavior, with ES resistivity seen at very small temper- 
ature and Mott resistivity at larger temperatures. This 
VRH is likely the result of some other source of disorder, 
unrelated to the positions of cations, which produces fi- 
nite DOGS near the Fermi level even dX v — 2. For 
example, if the NC diameters are not uniform, but are 
drawn from some distribution with finite width, then the 
energy levels corresponding to the IS and IP states are 
smeared. If the distribution of NC diameters has wide 
tails, then the IS and IP energy levels can be smeared 
as far as the Fermi level, producing a finite DOGS near 
the Fermi level, as shown schematically in Fig. [TUl 

The overlap between some IS and IP energy levels 
produces rare NCs with n = 3 or n = 1 whose energy is 
very close to the Fermi level. Such rare, mobile electrons 
are free to rearrange themselves in order to satisfy the ES 
stability criteria, and in doing so they produce a small 
Coulomb gap at the Fermi level (see Fig.lTU]). As a result, 
the resistivity follows the ES law at very small T, and 
the Mott law at larger T, where the DOGS sampled by 
electron hops is essentially constant. This is precisely 
what is seen in experiment [7|. 

It is worth mentioning that ionic liquid gating of NC 
arrays allows one to measure the total electronic charge Q 
as a function of applied gate voltage, or, in other words, 
the differential capacitance of the array C = dQ/dV. In 
arrays of small spherical NCs, where the quantum gaps 
A dominate over Coulomb energies, most electrons en- 
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FIG. 10. (Color online) Schematic picture of the density of 
states at i^ = 2 in the presence of fluctuations in the NC di- 
ameter D. a) If D has some wide-tailed distribution, then the 
IS and IP energy levels are broadened and have a finite over- 
lap, b) Spatial correlations between rare IS and IP energy 
states near the Fermi level produce a Coulomb gap, so that 
ES resistivity is seen at very small temperatures and Mott 
resistivity is seen at larger temperatures. 



VII. CONCLUSIONS 



In this paper we have used a simple theoretical model 
and a computer simulation to show how both activated 
transport and VRH arise in arrays of doped semicon- 
ductor NCs. Our primary result is illustrated in the 
phase diagram of Fig. [2j when the doping level v and 
the quantum confinement energy A are sufficiently large, 
and when the temperature T* is sufficiently small, the 
resistivity of the array is characterized by ES VRH. Such 
VRH is driven by the fluctuations in donor number from 
one NC to another, which lead to spontaneous charging 
of NCs as electrons depopulate higher quantum energy 
shells and fill lower ones. 

We have also identified a striking feature of the DOCS 
in NC arrays: the presence of "reflected Coulomb gaps" 
at electron energies ±2e^/KZ?, which are a consequence of 
the ES stability criteria and the discrete charging spec- 
trum of NCs (see Fig. |4]). This feature is even more 
prominent in large NCs with external impurity charges 
(Fig.®. 

The effect of additional disorder, such as fluctuations 
in NC size, remains yet to be explored quantitatively. We 
conjecture, however, that for chemically doped NCs our 
results will be largely unaltered by the addition of such 
disorder. For the case of NCs gated by ionic liquid, this 
external disorder seems crucial only for explaining the 
presence of Mott VRH at particular values of v (see Fig. 
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In our recent paper |28l|. we studied the hypothetical case 
where cations are large enough that only one cation can 
enter a pore in the NC array. In this case, due to the 
Coulomb interaction, the cations form a crystal struc- 
ture within the pores of the crystalline NC array. This 
situation is different from the model considered in Sec. 
rvl where ions are small and are introduced at relatively 
large temperature. In Ref. 12^ we argued that in the for- 
mer case the peak in capacitance corresponding to the IS 
shell splits into two delta-function-like peaks, such that 
one electron enters every NC at two particular values of 
the voltage. 



